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Abstract 

We develop a complete and consistent description for the hadron spectra from 
heavy ion collisions in terms of a few collective variables, in particular temper- 
ature, longitudinal and transverse flow. To achieve a meaningful comparison 
with presently available data, we also include the resonance decays into our 
picture. To disentangle the influences of transverse flow and resonance decays 
in the mr-spectra, we analyse in detail the shape of the my-spectra. 

1 Introduction 

Our current understanding of QCD results basically from high energy experiments 
with small collision systems suffering hard interactions which are relatively easy to 
analyse. In a first attempt to test QCD predictions for larger systems, especially 
the predicted phase transition from hadronic matter to the hypothetical quark gluon 
plasma, existing accelerators were modified to experiment with nuclei instead of only 
protons. The first round of experiments with nuclear beams took place during the 
years 1986-1990 at the AGS (Alternating Gradient Synchrotron) of the Brookhaven 
National Laboratory (BNL) and at the SPS (Super Proton Synchrotron) at CERN. 
While the biggest possible projectile nuclei 2S Si (BNL) and 32 5 (CERN) do not 
deserve the title "heavy-ions", the situation on the experimental side will consider- 
ably improve with the Au-be&m at the BNL, which is already in operation, and the 
planned P6~beam for CERN which is planned for 1994. To fully utilize the new pos- 
sibilities new methods for analysing the heavy-ion data have to be developed, which 



1 



are capable to characterize the physical situation over the whole range of nuclei and 
for both the energies at BNL and CERN. 

In this paper we want to develop a phenomenological model for the hadronic 
matter by starting out with thermalization as the basic assumption and adding 
more features as they are dictated by the analysis of the measured hadronic spectra. 
Eventually we will show that under the assumption of collective flow almost all 
hadronic spectra can be described by the model. The model can be extended to a 
consistent hydro dynamical description Q, which provides the link to the possible 
initial conditions and can in turn be used to extract new information about the final 
state of the hadronic system g, [| . 

We will start by explaining our choice of data in section ||. In section |3] we 
define the stationary thermal model for particle emission and show its failure for the 
measured ir~ m-p-spectra. We subsequently refine it by introducing resonances and 
their decay contributions to the m/r-spectra of all measured hadronic species which 
can then be described successfully by a uniform temperature. The analysis of the 
rapidity distributions in section 4.1 leads us to the introduction of longitudinal flow. 



On the other hand, as shown in section 4.2, the existence of transverse flow can not 



be established from the data. However this can be clarified by a closer theoretical 
investigation of the reaction prior to freeze-out, which was shortly reported in M 
and will be presented in detail in 0. Finally in section |E| we critically assess the 
relevance of our model in the light of pp data and give a conclusion of our work. 



2 Choice Of Data 

From the large amount of data from many experimental groups at the BNL and at 
CERN, which have been measured with many different targets and projectiles from 
27 Al to 184 W Q, we want to pick out one set of data which is for our purpose easiest 
to interprete: Transverse mass and rapidity distributions of pions ||, protons M, 
Kg, A und A Q, as measured in central collisions of 32 S with S at 200 GeV/nucleon 
by the NA35 Collaboration at CERN with a streamer chamber. 
For our selection we have the following reasons: 

• The high beam energy at CERN separates kinematically the central reaction 
zone around rapidity y cm = 3 from the fragmentation regions of the target and 
projectile at y = and y = 6. At BNL the span between target and projectile 
is smaller (y pro j — Viab = 3.4), and the width of each of the fragmentation 
regions Ay ~ 1 — 2 leads to a mixing of all zones. 

• The higher energies at CERN also generate a stronger longitudinal expansion 
of the matter, requiring at least cylindrical symmetry for the model, which 
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we will introduce in the next section, to describe the data. The BNL data 
for the pions Q] are closer to an isotropic momentum distribution and might 
be described well by a spherically symmetric model as we have presented one 
already some time ago || . 

• The streamer chamber from NA35 provides for a large number of measure- 
ments in the same experiment and thus guarantees easy comparibility. In this 
paper we will analyse the spectra of pions, protons, K®, A and A, which form 
the bulk of the hadronic matter. 

• The symmetric collision system S+S together with the selection of the events 
with 2% of the highest multiplicities minimizes the number of spectator nucle- 
ons, which did not collide with other nucleons at all. We are not interested in 
these nucleons, because they are similar to the majority of the nucleons in pA 
collisions, which suffer only a small energy or momentum transfer each. We 
are mainly interested in the participant nucleons, which generate the highest 
energy densities and where we can expect the strongest collective effects and 
would look for a quark gluon plasma. 

• The kinematic symmetry is experimentally fortunate, because only the wider 
open, easier accessible rear half of phase space (y < 3) has to be measured, 
the other half can be obtained by reflection around y = 3. E.g. by measuring 
the negative hadrons in the interval 0.8 < y < 3.0 and reflecting around y = 3, 
NA35 covers almost the whole rapidity gap between projectile and target. 
Also the isospin symmetry of 32 S can be employed for the equality of proton 
and neutron spectra, as well as 7r _ , 7r° and 7r + , allowing for a determination 
of the proton spectra by subtracting the spectra of negative tracks from the 
spectra of the positive tracks with only a small error introduced by the possible 
inequality of the K + and K~ spectra and p contributions ||. 

These advantages outweigh the disadvantage that S+S is not the biggest possible 
collision system to date. It contains with full overlap at most 64 nucleons, whereas 
a bigger target, e.g. 197 Au, would increase the number of participating nucleons to 
115 in a simple geometric picture and moreover would raise the baryon and energy 
densities as well. However, besides the blurring of the theoretical explanation by 
the big number of non-participating nucleons, for an asymmetric collision system the 
phase space has also to be measured in forward direction, which is experimentally 
increasingly difficult for the heavier collision systems. These data are thus not 
available yet, 32 S+^ W7 / im Ag is in the analysis right now M . The biggest asymmetric 
system with a complete set of published hadron data is up to now 16 O + 197 Au with 
a maximum number of participants (approx. 70 from geometry) again similar to 
S+S. 
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Since the fundamental dynamics of an ultrarelativistic heavy-ion reaction will 
basically be determined by strong interactions, the regular hadrons (pions, protons 
and kaons) carry the major part of the total energy jn]]. These particles spectra 
are most important for the global picture of the reaction. To a large extent we will 
rely on the pion spectra, because the pions are by number and also by energy (2/3 
of total) the most important group and can be measured best, i.e. NA35 determines 
the 7T~ spectra by recording all negative tracks and only correcting them for electron 
tracks. The contamination by other particle species, i.e. by antiprotons or K~ , is 
quantitatively small and can be accounted for on a statistical basis. 

In this paper we will only be interested in the shape of the spectra. The absolute 
multiplicities or the particle ratios convey little direct information about the reaction 
dynamics and require additional knowledge about the volume of the interaction zone 
or the equation of state, respectively. We will defer these questions to the more 
detailed hydrodynamical simulation of the collision zone in . 

3 Spectra From A Stationary Thermal Source 
3.1 A Purely Thermal Source 

In order to clarify the notation we start with the invariant momentum spectrum of 
particles radiated by a thermal source with temperature T: 

E fn = dn = _gV^ E ^e^/t (1) 

d 3 p dymTdrriT d<p (27r) 3 

where g is the spin/isospin-degeneracy factor for the particle species and fi the grand 
canonical potential [i = bfif, + s/x s as originating from its baryon and strangeness 
quantum numbers b and s. For simplicity we neglect quantum statistics with the rea- 
soning that its influence will be rather small at the low densities where the particles 
typically decouple from each other and where the spectra are computed. This also 
holds true for the pions, whose Bose-character is only then important for the spec- 
tral shape when chemical potentials close to the pion mass are introduced |jll| via 
non-equilibrium arguments, which are somewhat beyond the scope of our analysis. 
V is the volume of the source, giving together with the factor e^^ T the normalization 
of the spectrum, which we will always adjust for a best fit to the data, because we 
are only interested in the shape of the spectra to reveal the dynamics of the collision 
zone at freeze-out. In the remainder of the text we will always give the spectra in 



terms of rapidity y = tanh -1 (jpi/E) and transverse mass my = ym 2 +Pj>- Also 
h = c = ks = 1. 

We obtain the transverse mass spectrum dn / '(mrdmr) by integrating over ra- 
pidity using the modified Bessel function Ki, which behaves asymptotically like a 
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decreasing exponential: 

dn \ ' 



V is ( m T\ m T ^T , m-r/T m\ 



mTdrriT 2-7T 2 \ T ' 




Fig. 1: Purely thermal 7r 
rriT-spectrum in comparison 
with data from NA35 S+S 200 
A GeV @. A fit of a purely 
thermal source (eq. ||) was at- 
tempted to the data with the 
temperature T as a free param- 
eter. The absolute normaliza- 
tion was adjusted best possible, 
nevertheless the deviations at 
medium and high nix are sig- 
nificant. 



While the basic characteristics of the measured m/p-spectra is indeed the expo- 
nential decay over several orders of magnitude with an almost uniform slope 1/T 
when plotted over itit, the ir~ spectrum from NA35 shows a significant concave 
curvature (Fig. [l]) on top of the exponential dropoff, which is not visible in the rriT- 
spectra of the other particles (p, K®, A, A). This curvature makes it impossible for 
a single thermal source to fit both the low and the high region at the same time. 
(It does not stem from the curvature of the K\ function, since that one is much 
smaller and in the opposite direction.) The effect has been known under the name 
"low-rwr enhancement" since the first data from the ultrarelativistic heavy-ion ex- 
periments became available §, and has stimulated many theoretical interpretations 
(for an overview see e.g. fl2|j). 

Before we will devote ourselves to that problem in detail we finish the discussion 
of the thermal source with the rapidity distribution, which results from integrating 
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the invariant momentum spectrum over the transverse components: 



dn-th V of m 2 m 2 2 \ / m \ 

~^ = Wf fe + Tc^ + ^h^) eXP l"T C ° Shy J • (3) 

which reduces for light particles to dn/dy oc cosh _2 (y — j/q) up to third order in m/T. 
This rephrases the fact that the rapidity distribution of massless particles from 
an isotropic source is always the same regardless of their momentum dependence 
(e.g. J|). Its full width at half height T^ hm ss 1.76 is in sharp contrast to the 
experimental value for the pions of = 3.3 ± 0.1 § (see dotted line in Fig. @). 

The particle mass further narrows the rapidity distributions, leading for m > T 
to an additional gaussian with the width parameter a 2 = m/T, and increases the 
deviations from the measured distributions. Clearly also here are improvements 
necessary, which we will attempt in section [|. 

3.2 Resonance Decays 

The attitude of the previous section was a bit too naive: Apart from the directly 
emitted pions of assumedly thermal origin the detector also detects pions, which 
originate from the decay of resonances, which are also generated in the reaction, e.g. 

p° — ► 7r + 7r~ , uj — > ir + ir°ir~ or A —> Nir~ 

and which have a very different spectral shape, thus distorting the straight expo- 
nential dropoff. 

The abundant production of these resonances is an experimental fact, which 



has been investigated closely in pp collisions [13], where because of the much lower 
multiplicities the individual resonances can be reconstructed. Using these data as 
a guidance for our implementation in nucleus-nucleus collisions, we find that the 
absolute multiplicities of the resonances are distributed like oc exp(— uir/T) with 
T ~ 180 MeV, justifying a thermal descripton as an approximation as long as the 
temperature is close to this value. For the baryonic resonances (A, A*, . . .) we also 
fix the baryon chemical potential to be fiB = 200 MeV and distribute the strange 
resonances (K*, A*, . . .) according to strangeness neutrality fll4j| . The spectral shape 
for the resonances has also been shown in pp experiments to be exponential with 
T ~ 180 MeV and can for the heavy-ion collisions quite safely be assumed to be 
thermal since frequent elastic collisions will equilibrate the resonances in the same 
way as the ground state particles. 

The decays of the heavier resonances into their lighter daughter particles has 
been investigated in detail by Sollfrank et al. in and others [16|. Since we will 
analyse the mechanism in detail later in App. |b], we will here only sketch the method 
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of computation. The spectrum of the daughter particles is obtained by integrating 
the momentum distribution of the resonances with a factor f(W,pR,p) describing 



the decay phase space [15| within the kinematic limits of invariant mass W, rapidity 



wM 2 v R +) 



y and transverse mass mx 

. " = [ dW 2 [ dy R { dm 2 TR f(W, PR ,p) , Hr . .(4) 

ay mTdrriT J J J dyn mTRdrriTR 

y R TR 

To compute the ir~ spectrum we evaluate (|j) for the two-body decays p, K~* , K°* , 
A, E° , A* as well as for the three-body decays of u and r\ numerically. 

With these assumptions the transverse mass spectrum of the pions is again only 
a function of the temperature T and baryon chemical potential //& and one can try 
a fit to the data. As shown in Fig. |2[ we now succeed in reproducing the spectrum 
over the whole range in tut, with the temperature T corresponding to the slope at 
high niT- The kinematics of the resonance decays result in very steeply dropping 
daughter pion spectra and raise considerably the total pion yield at low tut- For the 
heavier p, A and A (Fig. ||) the change in the spectra is less pronounced, since 
the mass difference between the resonance and the ground state is much smaller. 
Quite remarkably, we realize that all independent fits to these five particle spectra 
seem to agree in T ~ 200 MeV within their statistical uncertainties Jul, 17]. 



We conclude that the resonance decays, as a necessary ingredient for any thermal 
model, provide already a good description of all measured hadronic my-spectra. The 
magnitude of the effect (2/3 of the pions come from resonances) as it was assumed 
here, is caused by the high fit temperature and has not been confirmed by further 
data, e.g. on the relative multiplicities of p, uj, A, etc. in S+S, so that there is 
still some room for other explanations (which should, however, always also take the 
resonances into account, albeit perhaps at a somewhat reduced level). 

Especially the possibility of a chemical potential for the pions in size of fi n ~ 
118 MeV at T = 164 MeV has been discussed [jlT| , which might originate from non- 
equilibrium effects and which would accumulate pions at low momenta because of 
the nearby divergence of the Bose distribution. With our treatment of the n~- 
spectra we cannot exclude this effect, but we can reduce its importance by first 
subtracting the always present resonance contributions. In a quantitative analysis 
of the pion multiplicity in our hydrodynamic model we determine that the non- 
equilibrium population of pions and the accompanying pion potential is much smaller 
than the pion mass and cannot influence the shape of the m^-spectra via the upwards 
curvature of the Bose distribution. 

For a further overview of the low m^-enhancement we would like to refer to |l2j . 
From our side the possibility of transverse flow was brought into the discussion with 
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Fig. 2: Thermal tt~ mr-spectrum with resonance decays in com- 
parison to data from NA35 S+S 200AGeV [|. The purely thermal n" and 
the two- and three-body decays of mesons and baryons add up to the total 
spectrum. The fit of the total 7r~-spectrum from a thermal source with the 
temperature T as a free parameter is very good over the full range. The abso- 
lute normalization have been adjusted for a best fit, /if, = 200 MeV has been 
fixed. These results from [fl5| are repeated here for later comparison within a 
coherent presentation. 
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Fig. 3: Thermal mx-spectra of p, K°, A and A including resonance 
decays in comparison with data from NA35 S+S 200 GeV/n @, |). The 
purely thermal spectra (dotted) and the resonance decays (dashed) add up to 
the total spectra (solid). The respective fits to the my-spectra are all very good 
and give similar temperatures. The absolute normalizations have been adjusted 
for a best fit, /it, = 200 MeV has been fixed. These results from from Jljj] are 
repeated here for later comparison within a coherent presentation. 
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a model based on spherical symmetry ||. Here, however, we will adopt cylindrical 
symmetry and analyse the longitudinal and transverse flows separately. While the 
longitudinal flow can be extracted from the data, the same uniqueness cannot be 
achieved for the transverse flow, though consistency arguments suggest its existence. 



4 Spectra From Flowing Sources 
4.1 Longitudinal Flow 

We have already mentioned that the rapidity distribution of 7r~ (actually nega- 
tively charged hadrons) as measured by NA35 is much wider than a stationary 
thermal source could possibly predict. Also the inclusion of resonance decays does 
not improve this picture, since in the same way as they lead to a steepening of the 
mr-spectra at low rriT they cause a slight narrowing of the y-distribution. 

Obviously the momentum distribution of the measured pions is not isotropic, it 
rather has imprinted on it the direction of the colliding nuclei. The boost-invariant 



longitudinal expansion model, as it has been postulated by Bjorken [18] can explain 
such an anisotropy already at the level of particle production, which subsequently, 
after sufficient rescattering, leads to a boost-invariant longitudinal flow of matter 
with locally thermalized distributions. The observed particle spectrum then results 
from the summation of the spectra of individual thermal sources which are uni- 
formly distributed in flow angle n. However, since the model has been formulated 
for asymptotically high energies, where the rapidity distribution of the produced 
particles establishes a plateau at midrapidity, we cannot apply it directly to our 
current energies, where already half of the total rapidity gap of 6 units is eaten 
up by the target and projectile fragmentation regions and consequently the pion 
distribution rather looks like a gaussian. 

We modify the boost invariant scenario to account for the limited available beam 
energy by restricting the boost angle rj to the interval (?? m i n , r/ max ) p9| . The rapidity 
distribution is then the integral over the uniformly distributed thermal sources (|3|) 
boosted individually by r/: 

dn f dn t h 



— (y) = J dri—j—{y-ri) (5) 



mm 



linn 



The transverse mass spectrum is not affected by this operation. 

We can now try a fit of (||) to the measured ir~ y-spectrum using r/ max = —n s 
as the single free parameter, because the distributions are not very sensitive to T. 
The rapidity distribution of negative hadrons from NA35 is measured over a large 
acceptance region (0 < px < 2GeV/c) and can be identified with the n~ , since the 
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Fig. 4: tt~ y-spectrum with longitudinal flow in comparison with 
data from NA35 S+S 200 AGeV 0, which have been reflected around y = 3. 
The maximal fluid rapidity 77 ma x = 1.7 has been obtained from a fit of eq. (||) 
to the measured spectrum. The absolute normalization has been adjusted 
for a best fit, the agreement with the data is very good. The temperature 
T = 220 MeV has been taken from the fit to the mr-spectrum and has only 
a negligible influence on the y-spectrum, similar to the inclusion of resonance 
decays (dashed line). For comparison also the purely thermal spectrum without 
any longitudinal flow is shown as a dotted curve. 
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Fig. 5: y-spectra of p, K®, A and A with longitudinal flow in com- 
parison with data from NA35 S+S 200AGeV §, §. The flow ?7 max = 1.7 
seen in the pions (Fig. ||) describes all produced spectra satisfactorily, only the 
protons still carry a large fraction of their initial collision energy. The abso- 
lute normalization was always adjusted for a best fit, the temperature and the 
resonance decays (dashed) have no influence on the spectra. For comparison 
also the purely thermal spectra without longitudinal flow are shown as dotted 



curves. 
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additional contributions from K~ andp are small (5-10%) and probably distributed 
more or less homogeneously over the whole region ||. The data points in the not 
accessible region of the streamer chamber above y > 3.4 can be substituted using 
symmetry with respect to y = 3. 

We see in Fig. || that 7? max = 1-7 fits the measured 7r~-spectra quite nicely. 
Again, the inclusion of the resonance decays does not qualitatively change the ir~ 
distribution since the resonances are uniformly distributed over the boost rapidity rj 
and have thus the same minor effect on the integral @ as they have on the individual 
thermal spectra (|3|). This would be different for inhomogeneities in the resonance 
distributions: for example a high concentration of A's in the fragmentation regions 
would manifest itself by additional pions from that region, i.e. bumps in the rapidity 
distribution [17|. 



Looking at the other hadronic distributions (Fig. |5|), we see our interpretation 
basically confirmed. However, the experimental proton spectrum is much broader 
than the computed one and has a dip at central rapidity; the measured protons 
obviously carry still a big amount (50%) of their initial collision energy [10], and 
the spectrum contains many protons near the target and projectile fragmentation 
regions, which have not suffered sufficiently many collisions to have become part 
of the central fireball. Because there is no clearcut separation between the central 
region and the fragmentation regions in the experimental y-spectrum, our model, 
which is crafted for the central region only, is bound to fail for the fragmentation 
zone protons. 

The produced strange particles are much better reproduced: the theoretical K® 
spectrum is only slightly too broad, the one for A's slightly too narrow and the one for 
A's still satisfactory with the exception of the central data point which also has the 
largest statistical uncertainty. The strange particles are good indicators for collective 
flow, because all are produced particles which were absent in the original nuclei and 
are thus not contaminated by a cold spectator component at leading rapidities as 
the protons are. Also because of their bigger mass their flow component is more 
accentuated against the random thermal motion as in case of the pions. Moreover, 
the disappearance of the characteristic polarization of the A Q hints towards at 
least one other collision in the medium, which justifies viewing them as part of a 
thermalized system. For the A the situation is more complicated and the success 
of our model might be only accidental. We can imagine that they are produced 
centrally in particularly hard individual nucleon-nucleon collisions and have not yet 
approached thermal and chemical equilibrium between annihilation and production. 
We could thus expect to see in the A distribution besides the flow component the 
imprint of the production process at central rapidities. However, a recent chemical 
analysis of the various particle ratios |20|l, which views the A's as an integral part of a 
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Fig. 6: The longitudinal flow 77 max 

can be extracted from the rapidity dis- 
tributions of the hadrons. With excep- 
tion of the protons, which carry still a 
large fraction of their initial beam mo- 
tion, a fit with 77 max as a free parameter 
gives for all produced particles indepen- 
dent of the assumed temperature a lon- 
gitudinal flow similar to the pions with 
?7max « 1.7 at T ps 200 MeV. 



thermally and chemically equilibrated fireball appears to be phenomenologically very 
successful and in favor of a thermalized picture which has already lost its memory 
of the production process. 

We can also try a best fit of the hadronic y-spectra individually by adjusting ry max 
for each particle species separately (Fig. ^). Again we arrive at the same value of 
Vraa,x = 1.7 ±0.3 independent of the temperature. Unfortunately, the ^-distributions 
of the heavier particles have less statistics than the pion spectrum and do not extend 
to as small rapidities; hence the estimates of the widths relies heavily on the last 
data points. Therefore we will prefer to continue our reasoning with the better 
determined pion value. 



4.2 Transverse Flow 



The stepwise extension of the initial stationary thermal model by including reso- 
nance decays and a longitudinal flow component has provided a phenomenological 
description of almost all spectra - dn/dy as well as dn/dmj> - of pions, protons, K® , 
A and A as measured in the experiment NA35. The few exceptions can be explained, 
they show the limits of our model. However, this seemingly complete picture of the 
dynamics of the hadronic matter has some serious intrinsic theoretical problems: 



1. On one hand at temperatures around 200 MeV the mean free path of pions in 
hadronic matter is much less than 1 fm, as we can estimate from the thermal 
distributions and averaged cross section of pions with nucleons and with each 
other [15|. On the other hand the size and lifetime of the reaction region is 
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several fm, as can be guessed from the radius of the 32 S nucleus and causal- 
ity. The length of the collision zone will be even larger because the large time 
dilatation factors of the relativistic longitudinal expansion stretch the system 
quite long already during the formation time 77 ~ lfm/c. Consequently, the 
pions cannot leave the interaction zone at T ~ 200 MeV without further colli- 
sions, the reaction region cannot decouple thermally and should by continuing 
expansion force the pions to cool down further. 

2. The longitudinal expansion in the boost-invariant scenario corresponds to 
a solution of the hydrodynamical equations in (l+l)-dimensions. In fact 
the 1-dimensional expansion presumably dominates initially because of the 
anisotropic initial conditions. But it is inconsistent to assume that a ther- 
malized system expands collectively in longitudinal direction without generat- 
ing also transverse flow from the high pressures in the hydrodynamic system. 
This should be computed by hydrodynamics in at least (2+1) dimensions, if 
azimuthal symmetry is maintained, otherwise (3+l)-dimensions, because the 
transverse expansion will increase the cooling rate. 

3. Already transverse velocities of about 0.3-0.6 c, which are moderate compared 
to the longitudinal expansion with j3i = tanhr/ max > 0.9 c, could change the 
transverse mass spectra considerably by enhancing the particle yields at high 
niT- 

We want to address these problems with a model which contains resonance decays 
and both longitudinal and transverse flow in order to put the discussion of the 
longitudinal and transverse spectra on a common basis. In this paper we will analyse 
how much information can be extracted from the data with such a model. Using our 
phenomenological analysis presented here we will extend in the next step [Q, [| the 
model towards a global hydrodynamical description which enforces the mentioned 
theoretical consistency requirements. We observe (similar to [[17]]) that agreement 
with all the presented data can again be achieved and we prove the existence of 
transverse flow by theoretical means || |3| . 

Not only theoretical necessities suggest the inclusion of transverse flow. There 
are further arguments for the actual existence of the phenomenon: 

1. At the lower BEVALAC energies (E/A < 2GeV) a sidewards flow of nuclear 
matter is observed [M] . It is generated by the squeeze-out of nuclear matter 



in the collision and has been predicted by hydrodynamic computations [22]. 
Though this directed flow is of different nature as our collective expansion flow, 
it suggests the relevance of hydrodynamics also at higher energies. 
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Fig. 7: tv~ m/p-spectrum with resonance decays and transverse flow 

in comparison with data from NA35 S+S 200AGeV ||. The total ir~ spec- 
trum consists of the original pions and the decay pions. Based on a parabolic 
transverse velocity profile with (3 S — 0.5 c the fit with the temperature T as the 
free parameter is quite good over the whole mr-range. The absolute normal- 
ization has been adjusted for a best fit, ^ — 200 MeV has been fixed to allow 
for a moderate baryon population. 
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Fig. 8: The m^-spectra of p, K°, A and A with decays and transverse 
flow in comparison with the data from NA35 S+S at 200 A GeV @ § . Based 
on a parabolic transverse velocity profile with [3 S — 0.5 c the individual fits with 
the temperature as the free parameter show good agreement with the data at 
temperatures similar to Fig. ^. The absolute normalization has been adjusted 
for a best fit, /i^ = 200 MeV has been fixed as in Fig. |?]. 
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2. In a simple picture the transverse flow flattens, in the region p± < m, the 
transverse mass spectra of the heavier particles more than for the lighter par- 



ticles [23, g]. The AGS data [24] show clearly this tendency for pions, kaons 



and protons. 

3. Event generators for nucleus-nucleus collisions show a statistical effect which 
can be interpreted as transverse flow. The statistical averaging of the veloc- 
ities of many particles in a small volume can exhibit the transverse flow at 
AGS energies [^f|. The VENUS event generator for nucleus-nucleus collisions 
shows also for CERN-energies an increase of the spectra at high tut, once the 
produced particles are allowed to rescatter |2fJ]. The same effect is produced 
in our model by the transverse flow. 

On the other hand we have to admit that from phenomenology alone there is 
no need for an additional transverse flow effect, since all the data have already 
been described nicely. Currently we thus will not be able to prove the existence 
of transverse flow from the data alone, only an indirect proof based on the phe- 
nomenological analysis together with hydrodynamical calculations with a consistent 
theoretical treatment of the freeze-out process can be given Q. But in the near 
future bigger collision systems can be expected with the ^4u-beam at the AGS and 
with the installation of the P6-injector at CERN, which will lead to an increase of 
the collective effects. If we are currently not able to prove transverse flow directly 
from the data, the chances will be improving tomorrow, and we should develop our 
tools in the meantime. 

We compute the spectrum by boosting the thermal sources now both in longi- 
tudinal and transverse direction. We describe the transverse velocity distribution 
Pr(r) in the region < r < R by a self-similar profile, which is parametrized by the 
surface velocity j3 s 

^ (r)=/3s G?) n • (6) 

With n we can vary the form of the profile. We choose customarily n = 2, because 
the quadratic profile resembles the solutions of hydrodynamics closest [0. Anyway, 
the form of the profile is not important for the analysis, as we checked with the case 
n = 1. Leaving the details of the computation to App. [A|, the resulting spectrum is 
a superposition of the individual thermal components, each boosted with the boost 
angle p = tanh -1 (3 r : 



dn 

oc 



.. l\ irrtwIo l^^) Kl l^ptl) (7) 
rriT arriT Jo v 1 ' v 1 ' 

From Figs. |?| and ^ we see that also with a moderate transverse flow (/3 S = 0.5 c, 
i.e. (P r ) = 0.25 c) very good fits to all hadron spectra can be obtained, again 
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treating the temperature as a free parameter. The temperature of the local sources 
as extracted from the data is again suprisingly uniform but now much smaller than 
without transverse flow (Figs. ^ and ||), thereby reducing the contribution from the 
resonances. A quantitative explanation of this effect will be given later in App. [A|, 
qualitatively the effect is similar to a blue shift as caused by a rapidly approaching 
source, shifting particles to higher momenta. The heavier the particles, the more 
they profit from the flow velocity. 
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Fig. 9: All fitpairs (T,/3 S ) 
of the 7T _ , K°, p, A and A 
mx-spectra. Every point 
on these curves gives a good 
agreement of the computed rap 
spectrum (with resonances and 
transverse flow) with the re- 
spective measured spectrum. 
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We condense the information from the my-spectra with respect to the transverse 
flow in Fig. || where we show all possible fit pairs (T,/3 S ) for all hadronic spectra. 
Without transverse flow {(3 S = 0) the intercepts of the curves with the axis show 
the fitted temperatures of stationary thermal emitters including resonance decays 
(Figs. H & |3|), where all temperatures cluster around the value T « 210 ± 20MeV. 
Surprisingly, despite their somewhat different shapes, the curves stay together in 
one band and only weakly concentrate in the region (3 S ~ 0.5 c. 

Intuitively one would rather expect the contrary: Because flow increases the 
particle energies according to their rest masses ((E) « {E)th + 



mo„2 



coll )i 



the heavier 

particles profit more from the collective motion than the lighter ones. Since both 
transverse flow and the local temperature determine the slope of the m^-spectrum, 
a bigger effect from flow for the heavier particles should result in faster dropping 
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Fig. 10: Quality of the (T,/3 S ) 
fits for the tt~ mx spectra. The 

contour lines designate the areas of 
X 2 /NDF = 1, 2, 4, 8, 16, and 32. fi b 
has been fixed at 200 MeV. 
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Fig. 11: Fit pairs of tt~ , K°, 
p, A and A rriT spectra omitting 
the decay contributions. In this 
unrealistic situation there would be an 
intersection around f3 s « 0.5 c, which 
would strongly favor a sizeable trans- 
verse flow. 



(T, /3 s )-curves from the fits for K®, p and A. The different shapes of the curves 
should make it possible to find one unique intersection point, representing T and 
(3 S of the collective system || . Unfortunately, this simple argument does not take 
the resonance decays into account, which modify especially the pion spectra strongly 
1 15]. We show in Fig. 11 the curves, which unrealistically omit the resonance decays, 
and from which we would erroneously deduce a transverse flow of (3 S ~ 0.5 c from 
the unique intersection of all curves. However, as we see in Fig. ^, reality is more 
complicated, the pion curve is considerably modified by resonance decays so that in 
the end its shape resembles the heavy particles' and all values < (5 S < 0.7 c seem 
to be possible within the experimental uncertainties. 

The A do not quite fit into the picture, they seem a bit colder than the other 
particles including the protons. An explanation in terms of an admixture of colder 
A from the fragmentation region is supported by the the rapidity distribution but is 
in some contradiction to the higher proton temperature (which however has exper- 
imentally the biggest error bars). In a similar manner one could justify the higher 
temperatures of the A which seem to originate more from the hotter center. 
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We can conclude that the data alone do not make a unique statement about the 
existence of transverse flow. Instead of a single point for (5 S there is a band with 
many pairs (T, (3 S ) . The phenomenological analysis has fully exploited the available 
data and has to stop here. Which point on these curves corresponds to the freeze-out 
point of a consistent hydro dynamical evolution will be studied in a seperate paper 
|2|, H| . The theoretical methods presented there can use the gathered information to 
back-extrapolate towards the early stages of a heavy-ion collision and to eventually 
fix the amount of transverse flow. 

5 Discussion And Conclusions 

Naturally there is some doubt whether collective effects are actually as important as 
it might seem from the analysis presented here, since the S nucleus is still very small. 
Comparing S+S to pp minimum bias data, which should resemble most closely the 
average nucleon-nucleon collision in S+S, we realize that the rapidity distributions 
are rather similar in shape. Accounting for the isospin symmetry, the width of the 
S'+S'-distribution turns out to be only about 10% smaller than from pp ||. |T(J . 

The transverse momentum spectra show bigger differences when plotted as the 
(normalized) ratio of SS/pp p7|| . Apparently there is an enhancement at low pt, 
which could result from the increased population of resonances, and a significant 
enhancement at high pt, which we would contribute to the bigger transverse flow. 
However, a high p^-enhancement has also been observed in pA collisions [^] (Cronin 
effect), although this effect sets in only above p± ~ 1.5 GeV. It has been interpreted 
as resulting from multiple collisions of the partons in the nucleus [^] and as such 
is a QCD-specific effect. We interpret this multiple scattering already in pA as an 
indicator that collective matter flow can be generated, since in the limit of many 
scatterings hydrodynamic behaviour will eventually result. 

It would be too impulsive to deduce from the apparent similarity of pp and S+S 
spectra that the collective interpretation is wrong, since a) pp is by no means an 
elementary collision system which we understand in sufficient detail to serve as the 
antipode of a collective system and b) only a few of the observed features of the 
spectra can be fully reproduced by these two radically different philosophies which 
share only a small set of common principles as local energy-momentum conservation, 
relativistic space-time picture, etc. For other observables, which are not directly tied 
to the dynamics, e.g. strangeness, we already see a big enhancement compared to 
pp, thus indicating fundamental differences between both collision systems. 

A careful systematic study of the ^-dependence from pp up to Pb+Pb should 
allow to decide this question by showing the increasing importance of the collective 
effects for the bigger systems relative to the individual scattering processes. Un- 
fortunately at the present time the only other symmetric collisions at high energies 
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are pp collisions |3(|, because, as we already mentioned, we cannot directly compare 
with asymmetric (pA) collisions. The data from the no longer operating ISR (Inter- 
secting Storage Rings) from light nuclei (dd, aa) (3l]] would also be very interesting, 
but they probably have to be reanalysed for a direct comparison. 

Nevertheless, from our experience with S+S we can already risk the prediction 
that the spectral shape from Pb+Pb collisions will show only gradual changes from 
S+S and that a careful analysis will be needed to uncover the interesting physics. 
How much influence the creation of a plasma would have on the spectra can not 
be answered in this context without additional theory, but we would already like to 
caution big expectations, since the plasma would be created in the early stages of 
the collision and any signal would be leveled by the following hadronic stages with 
their own dynamics. 

In conclusion we have shown that a thermal model is perfectly possible for S+S 
collisions despite (or because of?) the similarity of S+S and pp spectra. The 
data force us to include resonance decays and longitudinal flow while they make 
no decisive statement about the existence of transverse flow. But through further 
theoretical investigations the quantitative analysis presented here can be used to 
prove the existence of transverse flow and infer other informations about the collision 
zone g|. 
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A Slopes From Transverse Flow 

Both transverse flow and resonance decays have a strong impact on the tut spectra. 
In this and the following section we will show the actual computation of the spectra, 
and we will try to disentangle both influences by analysing the slope of the rriT- 
spectrum. 

We follow the spirit of the boost invariant scenario by first defining the transverse 
velocity field in the central slice (using the boost angle p = tanh -1 (3 r ) for azimuthal 
symmetry 

u' u (t, r, z = 0) = (cosh p, e r sinh p, 0) (8) 
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and later boosting it in longitudinal direction (boost angle rj) to generate the whole 
velocity field [32] 



u M (p, rj) = (cosh/? cosh 77, e r sinhp, coshp sinhr/) . (9) 

Note that u M is not symmetric with respect to the boost angles p and rj. This results 
from the fact that Lorentz transformations for different directions do not commute. 
Having defined the velocity field, the invariant momentum spectrum is given by 

m- 3 

E ^ = jj^P)P Xd ^ « 7A3 / e-^-^ T p x da x , (10) 

where f(x,p) is the invariant distribution function, which we assume to be an 
isotropic thermal distribution boosted by the local fluid velocity u^, and we approx- 
imate the respective Bose and Fermi distributions by the Boltzmann distribution. 

This mathematical formulation measures directly the particle flow through the 
given hypersurface a as if the virtual walls of the fluid cells have suddenly disap- 
peared and the particles are flying isotropically in all directions, a should not be 
visualized as a real surface which emits radiation only to the outside. Instead it 
defines the borderline between hydrodynamical behaviour and free-streaming par- 
ticles, which are both idealizations, as a mathematical construct. More intuitive 
approaches lead to formulas which no longer obey the conservation laws. In re- 
ality the freeze-out hypersurface might be defined by the points of the last interaction 
of each individual particle and thus acquires a thickness of the order of the mean 
free path. This seems to be difficult to implement consistently with the dynamical 
evolution and for our purposes the simple model seems to be quite sufficient. 

We parametrize the hypersurface a(r, (f>, £) in cylindrical coordinates < r < R, 
< 4> < 2-7T and — Z < £ < Z in longitudinal directions. We invoke instantaneous 
freeze-out in r-direction to keep matters simple, but allow for more general shapes 
(t((),z(()) in longitudinal direction (e.g. hyperbolas |l8|) because of the large lon- 
gitudinal time dilatation effects: 

^{r, <f>,0 = ( t(C), r cos <p, r sin 0, z{C) ) (11) 



p^da 



I-'- 



'dz dt 

E d(~ PL dC 



rdrd(f>d( . (12) 



For the argument of the Boltzmann distribution we need the momentum of the 
particle in the center-of-fireball system 

u^p^ = rriT cosh(y — rj) cosh p — pr sinh p cos((p — if) . (13) 
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Because of azimuthal symmetry we can integrate over <f> making use of the mod- 



ified Bessel function Iq (z 

d 3 n 



(2^)- 1 / ( 
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(2vr)2 
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rdr exp 
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rriT sinn y— 



(14) 



mr cosh p cosh (y — rj) — p\ ( pt sinh/? 



For the transverse mass spectrum we integrate with the help of another modified 
Bessel function K\(z) = f£° coshy e~ zcosh ^ y dy: 
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(15) 



We see that the transverse mass spectrum factorizes, as long as temperature 
and transverse flow are independent of the longitudinal position in a longitudinally 
comoving coordinate system. There is only a factor Zx affecting the normaliza- 
tion. The invariant momentum spectrum factorizes only for small transverse flows 
(cosh p ~ 1) into a longitudinal and a transverse part. The rapidity distribution is 
almost independent of transverse flows, even if the latter becomes large, as we have 
checked explicitly. 

The general formula for azimuthal symmetry allowing non-instantaneous freeze- 
out in r-direction (dt/dr ^ 0) is rather similar, involving also Kq and I\\ 
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(16) 



We can achieve a better understanding of the transverse mass spectrum in ( ]T5[ ) 
by analysing its slope. Fixing the variables r, r/, T at reasonable values we focus on 
the rriT dependent parts. The slope in a semilogarithmic plot is 
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(17) 
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Fig. 12: m,T-slopes with 
transverse flow for pions, 
kaons and protons. Every 
group of curves corresponds to 
a fixed transverse expansion ve- 
locity f3 r . The approxima- 
tions of the curves by the pure 
exponentials (horizontal lines) 
are remarkably good already 
at low rriT. For comparison 
we included also the slopes of 
thermal sources (in Boltzmann 
statistics) with correspondingly 
blue shifted temperatures (solid 
curves). Compare also with 
Fig. " 
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We can immediately derive the limiting case for large m-r, since there mx/pT 
1, Kq/K% — » 1 and for finite flow (sinhp > 0) also — ► 1: 



lim 



d 



■In 



dn 



cosh p — sinh p 




Pr 



(18) 



m T -*oo drriT ~ ^rriTdmT' T T y 1 + j3 r 

The apparent temperature, understood as the inverse slope at high m-r, is then 



larger than the original temperature by a blue shift factor, see also p5fl : 



f eff 



'l + A- 

l-Pr 



(19) 



For small values of rriT the situation is much less clear, as one can see from 
Fig. (12|. We cannot expect a simple expansion around rriT = oo to deliver fully 
satisfying results, because apart from the shape of the Bessel functions also the 
particle mass mo (via pt = y/rriT 2 — rriQ 2 ) and the fluid rapidity p influences the 
slope. Nevertheless, the basic tendency can be discussed with this expansion 
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(20) 
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The first term is again the blue shift factor, the additional terms are responsible for 
the curvature of the spectrum. The deviations from the exponential behavior are 
largest for large flows (sinhp S> 0) and have the biggest effect at low rriT, where the 
shape of the Bessel functions tends to steepen the spectra, while big particle masses 
strongly flatten the spectrum there ||. 

One can extract regions from the figure, where the slope approaches the limiting 
value satisfactorily. For pions and kaons we have m-r > 1000 MeV, whereas for 
nucleons in the case of strong flow the limiting slope will be approached only above 
rriT > 2000 MeV. This has to be contrasted with the behavior of the resonance 



decays which only affect the low my region (see Fig. 13). 

Since in a realistic picture there are various components of transverse flow, we 
account for them by a superposition of simple exponentials, each with its own blue 
shift, denoted by b r ((3 r ) = ^/(l — /? r )/(l + f3 r ), to compute the blue shift factor of 
the integrated spectrum. 

d . r R 

lint 



drriT 

The integrals could be solved analytically for n = 1 and 2, but it is more instruc- 
tive to expand them into a power series in 1 — b s under the assumption that ttit/T 
is not too big: 



1-^(1-6.) -^P^(l-M 2 + ... (n = 2) 
l-|(l-b s )- 2(m % T) - 5 (l-M 2 + ... (n=l) 

1 -(l-b s ) (n = 0) 



(22) 



where the terms in third order do not improve the approximation considerably. 
Approximating the blue shift factor to first order by b s ~ 1 — f3 s , we arrive at a 
surprising coincidence with the average fluid velocity 



V 1 + \Pr) n + 2 

which is a handy way to estimate the behaviour of the fits in Fig. ^. 

In conclusion we find from our simple quantitative estimates that the influence 
of the transverse flow on the transverse momentum spectrum is strong. We also 
investigated the influence of the geometry of the freeze-out hypersurface on both 
the y- and on the m^-spectra. General formulas developed in |36|, show that 
the shape of the hypersurfaces influences the spectra via the partial derivatives of a 
with respect to r and z. Similarly to the spherical case || we have not encountered 
big differences in the width of the y-spectra as well as the slope of the mj- spectra 
as long as we confined ourselves to realistic freeze-out hypersurfaces. 
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Fig. 13: my-slopes of the 
7T~ with resonance decays. 

For realistic values of T and 
Pb we show the slopes resulting 
from the addition of always one 
resonance channel to the purely 
thermal 7r~-spectrum. At high 
rriT only the p decay changes 
the slope by a few percent. At 
lower m,T many individual con- 
tributions (especially A and u) 
add up to a significant modi- 
fication of the total spectrum. 
Interesting is the comparison 
with Fig. [12] (scales are differ- 
ent). 



B Slopes From Resonance Decays 

Computing the spectrum of the resonance products is generally a complex task 
and requires a computer to do the phase space integrals numerically. However, 
investigating the simpler case of the two-body decay analytically can give us some 
insights into the resulting spectra. We start from the general formula in fTql , now 
specialized for 2-body decays: 



d 2 n 



dymTdrriT 



m R b 



y ( R +) 



m| cosh {y - Vr) - Pt 



(24) 
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By introducing a thermal spectrum for the resonances which is independent of 



27 



rapidity d m™dm T = Ci m TR,e mTi? / T , we can solve one further integral analytically: 



y { R ] 

d 3 n m R f dy R 



(25) 



dym T dm T 4vrp* ^ ^ cosh 2 (y _ ^ _ p 2 

yR 

xirde-o/r { 2 (a 2 + 6 2 )/ (A) _ (4a6 + 2Tb)h(±) } 

where a = \{rnf R + Jn^) and b = ^(ra^ — rrSp R ). For large mr a and b are 
approximately independent of y R , 

m R E*m T _ , m R p*p T z /oc n 
a « k := a , o ~ k — := b (lb) 

and we can extract in first order the slope of the transverse mass spectrum since 
the Bessel functions Ii(b/T) behave asymptotically like exponentials exp(6/T). The 
resulting exponential is exp((6 — d)/T) giving a slope at high rriT 

d ln dn _ m R (E*-p*) l ^ 



drnj- niTdrnT m 2 

E* and p* are energy and momentum of the daughter particle in the rest system 
of the resonance and are for the two body decay already given by the particle masses 



(p* is listed for every decay channel in The estimate (27) coincides very well 

with the slopes of the numerically computed spectra at large mr- 

Further approximations reveal for the slope ^»?>f or th effective temperature 



T eff « ^—T (28) 
m R 

Since p* < m R the spectra of the daughter particles are generally steeper than the 
original spectrum of the resonance, which translates into a lower apparent temper- 
ature of the pion spectra at low ttit (Fig. |l^) . 

We can also estimate the width of the rapidity distribution by looking at the 
integration limits for y R . The maximal rapidity difference between resonance and 
daughter particle is restricted by kinematics to 

\VR ~y\< sinh" 1 (^-] (29) 
\m T J 

For the A-decay p* = 227 MeV, so that the decay pion can move at most 1.27 
units of rapidity (at rriT = thq), and the nucleon at most 0.24. The actual (y R — 
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y)-distribution of the daughter particle will naturally be much smaller, as can be 
inferred from the width of an isotropic distribution, which is for massless particles 
r fwhm _ 2 x 0.88 (eq. g). This similarity to a thermal source also explains the small 
influence of the resonance decays on the shape of the rapidity distribution (Fig. Q), 
which is only scaled by the decay contributions. 

In conclusion we have shown why the decay spectra are concentrated in the low 
m,T region. However, because of the superposition of many channels and the complex 
structure of the decay spectrum itself, we cannot estimate the influence by an easy 
method and have to resort to the full numerical evaluation of (||). 
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